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ABSTRACT 

We present a simple analytic model for the structure of non-relativistic and relativistic radiation mediated 
shocks. At shock velocities /3j = Vj/c > 0.1 the shock transition region is far from thermal equilibrium, since 
the transition crossing time is too short for the production of a black-body photon density (by Bremsstrahlung 
emission). In this region, electrons and photons (and positrons) are in Compton (pair) equilibrium at temper- 
atures r, significantly exceeding the far downstream temperature, ^ « 2(en,fl^c^y^'^. > lOkeV is 
reached at shock velocities /3s ~ 0.2. At higher velocities, /3s > 0.6, the plasma is dominated in the transition 
region by e± pairs and 60 keV < < 200 keV. We argue that the spectrum emitted during the breaking out of 
supernova shocks from the stellar envelopes (or the surrounding winds) of Blue Super Giants and Wolf-Rayet 
stars, which reach f3s > 0.1 for reasonable stellar parameters, may include a hard component with photon en- 
ergies reaching tens or even hundreds of keV. Our breakout analysis is restricted to temperatures < 50 keV 
corresponding to photon energies hiy < 150 keV, where pair creation can be neglected. This may account for 
the X-ray outburst associated with SN2008D, and possibly for other SN-associated outbursts with spectra not 
extending beyond few 100 keV (e.g. XRF060218/SN2006aj). 

Subject headings: shock waves — radiation mechanisms: nonthermal — X-rays: bursts — supernovae: 
general : individual (SN 2008D) 



1. INTRODUCTION 

During core collapse supernova explosions, strong shock 
waves traverse the exploding stars' mantle/envelope. These 
shocks move with high velocities, /3j = Vj/c > 0.01, through 
dense, n ^ 10'^ to 10^" cm"^, optically thick plasma. Under 
such conditions, the shocked plasma's thermal energy density 
is dominated by radiation as long as 
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where obb = 7r^/15(/!c)"^ is the Stefan-Boltzman energy den- 
sity coefficient. In such shocks, the photons that are gener- 
ated in the downstream region diffuse upstream and decelerate 
the incoming plasma by colliding with the electrons, which 
in turn stop the nuclei through collective plasma interactions 
(e.g. an electrostatic field). The typical width of such Radia- 
tion Mediated Shocks (RMS) is the length photons can diffuse 
into the upstream Ij ^ I /nu/3s<JT, where ctj- is the Thompson 
cross section, and the subscript u denotes in this paper the 
upstream values of parameters. 

The structur e of non-relativ istic radiation medi ated shocks 
was studied bv lColgatd (119741) and lWeaveri (11976b . motivated 
by the idea that the post shock plasma can reach high tem- 
peratures, T > \Q MeV, leading to the generation of deu- 
terium, as an alternative to Big Bang nucleosynthes is. Non 
relativistic numerical calculations were preformed bvl Weaveil 
(119761) for shock velocities /3s < 0.3. Self-Consistent solu- 
tions of the hydrodynamic, radiation generation (dominated 
by Bremsstrahlung) and diffusion equations were obtained, 
confirming that RMS are consistent physical structures. At 
low velocities, /3s < 0.1c [c.f. equation (fl4l ll thermal equi- 
librium is maintained throughout the shock. In such con- 
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ditions the shock temperature monotonically increases be- 
tween the far upstream and the far downstream values (e.g. 
IZel'dovich & Raized 119661) . For higher velocities, it was 
found that the plasma in the shock transition region de- 
parts from thermal equilibrium, and the electron tempera- 
ture reaches values of tens of keV, greatly exceeding the 
far downstream values (but much smaller than the suggested 
10 MeV temperatures). The results of the non-relativistic cal- 
culations are not applicable to shocks velocities exceeding 
/3s ~ 0.2, where the high temperatures reached in the shock 
transition require inclusion of relativistic processes (e.g. pair 
production), and where the evolution of the photon distribu- 
tion can no longer be simply described as diffusion in space 
and momentum. 

Detailed studies of Compton scattering in the converg - 
ing flow of n on-relativistic RMS fBlandford & Pavnd (11981b . 
lLyubarskir& Sunyaev (1982), Becker (1988) and recently 
in the context of shock-breakouts, Wang et al. ( 2007)] have 
shown that photons traversing through RMS can significantly 
increase their energy while scattering back and forth in the 
converging flow in a manner similar to Fermi-acceleration of 
Cosmic rays. The energy photons can reach by such 'bulk' 
acceleration is Umited to 



hi^^,, ~ 0.2m,c'/3^ ^ KA/O.iy keV, 



(2) 



above which the energy loss per collision due to the Comp- 
ton recoil is larger than the energy gain due to the converg- 
ing bulk flow. As we show here, for velocities /3s ^0.1 the 
electrons (and positrons) are "heated" in the immediate down- 
stream (the region where the photons are generated) to ener- 
gies exceeding those given in Eq. ^ [c.f. Eq.(fT8]ll. Thus, 
photons are generated with typical energies larger than given 
in Eq. implying that 'bulk' photon acceleration is only 
relevant for lower shock velocities, /?, < 0. 1 , and for sub- keV 
photons (and for densities low enough so that the downstream 
thermal temperatures are considerably smaller than /it'max)- 
In this paper we present a simple analytic model for the 
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structure of radiation mediated shocks. The model a ccu- 
rately reproduces the numerical results of lWeaverl (Il976h for 
Ps ^ 0.2, and provides an approximate description of the 
shock structure at larger velocities, (3^ I. A detailed com- 
parison of the results of the simple model presented here 
with the results of exact solutions of the shock structure 
(self-consistently solving the photon transport equation and 
the plasma flow equations for mildly and hig hly relativistic 
shock s) will be presented in a followup paper jBudnik et al. I 
HoT^. We assume in our calculations that photons are gen- 
erated mainly by Bremsstrahlung emission. This assumption 
is valid in the absence of strong magnetic fields (note that di- 
rect bound-free transitions do not lead to net generation of 
photons while the contribution of bound-bound transitions is 
negligible under the conditions relevant to shock-breakout). 

We first present in § |2] an analytical model for Non- 
Relativistic RMS, deriving simple analytic approximations 
for the non-equilibrium temperature and for the width of the 
temperature profile, and then extend our model in §[3] to rel- 
ativistic RMS. In section § |4] we discuss the implications 
of our results to the expected properties of SN X-ray out- 
bursts. We show that breakout velocities > 0.1 may be 
reached in the explosions of B lue Super Giant (BSG) and 
Wolf-Rayet (WR) stars (§ |4!2] i, and argue that the X-ray 
outbursts accompanying breakouts from such stars may in- 
clude a hard component wi th ph oton energies reaching tens 
or even hundreds of keV (§ 14.31 ). The implications of our re- 
sults to the interpretation of recent X-ray outbursts associated 
with SNe (XRO080109/SN2008D, XRF060218/SN2006aj) 
are discussed in § 14.41 We summarize the results and discuss 
their implications in § |5] 

2. NON-RELATIVISTIC RADIATION MEDIATED 
SHOCKS 

Consider a radiation mediated shock in a plasma consist- 
ing of protons, electrons and photons. The far upstream con- 
sists of cold protons and electrons moving with velocity PsC 
in the shock rest frame, while the far downstream consists of 
protons, electrons and photons in thermal equilibrium with a 
temperature Tj. 

In this section we construct an analytical, approximate 
model for the steady state profile of such a shoc k. W e start 
by writing down the downstream conditions in § 12.11 Next, 
we discuss the velocity profile in § 12.21 and argue that the de- 
celeration scale is roughly (/35n„(T7-)~'. In § 12. 31 we show that 
for large shock velocities, (3^ Jc^ 0.1, the temperature can have 
a non trivial profile extendi ng to distances much larger than 
the deceleration length. In § 12. 41 we estimate the temperature 
in the shock transition and show that it can greatly exceed T^. 
We show that our simple analytic results are in good agree- 
ment with detailed numerical calculations. Finally, we use the 
results derived in § |2.H -§ |2!4l to describe in § |2.5l the structure 
of fast NR radiation mediated shocks. 

2.1. Downstream conditions 

The energy and momentum flux are dominated by the pro- 
tons' momentum and kinetic energy in the upstream and by 
the photons in the far downstream. Proton number, momen- 
tum and energy conservation require that 

P-iJ = nuPsiPs- Pd)mpC^ , 

^P^jPd = riMP] - fil)mpc'l2, (3) 



where the subscript d denotes the downstream values of the 
parameters and we used the fact that the energy density in 
photons is = 3py. Equations (O imply fid = Ps/'^ ■ Thermal 
equilibrium implies 



e^,d - o-bbTci ■ 
The far downstream temperature is thus 



(4) 



/315 \^''^ 
Td= f ^«„£7zVj «0.16(ei«,u5)'/^keV, 



.0.13«:/;^<^keV 
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where «„ = 10'^«„ ig, e = 0.5l3^mpc'^ = lOei MeV and /?, = 
O.lA.-i- 

2.2. Deceleration 

Next we estimate the proton deceleration length scale. 
Once a proton reaches a point in the shock where the energy 
is dominated by photons, it experiences an effective force 

A^^MpC^ ~ ^7-^,6^ ajnuP^MpC^, (6) 
implying a deceleration length of 



^dec — /^s 



d/3 
dx 



CTTnuPs 



(7) 



We note that once the photons dominate the pressure, they 
cannot drift with the protons, as this will imply an energy flux 
greater than the total energy flux. Thus, the drag estimated 
in equation (|6]l is unavoidable. The energy and momentum 
of the protons are thus transferred to the photons on a length 
scale Ldec iiTespective of the details of the mechanisms that 
generate the energy in the radiation field. We note that the 
transition length, Ldec, is of the order of the distance a pho- 
ton can diffuse upstream before being advected with the flow. 
The deceleration of the fl ow can be solved analyticall y (e. g. 
iBlandford & Pavnell98lL and references within, cf Eq. ( |A6l l). 
For the strong shocks considered here. 
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where x is the distance measured in the shock frame and the 
point X = corresponds to /3//3rf - 1 « 0.25. 

2.3. Temperature profile length scale 

The region of the shock profile over which the temperature 
changes before it reaches Td can be extended to distances that 
are much larger than Ljec- 

To see this, consider the length scale that is required to gen- 
erate the density of photons of energy ^ Td in the downstream, 
determined by thermal equilibrium, n-^.eq ~ P-^j/Td, 



' Pc 



/t^y eq 



(9) 



where 2-y_eff. is the effective generation rate of photons of en- 
ergy 3Td. We use here the term "effective generation rate" 
due to the following important point. Photons that are pro- 
duced by the emission mechanism at energies <C Td may still 
be counted as contributing to the production of photons at Td, 
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since they may be upscattered by inverse-Compton collisions 
with the hot electrons to energy ^ 7^ on a time scale shorter 
than that of the passage of the flow through the thermalization 
length, Lj / PdC- Note, that the relevant Compton y parameter. 



AT AT 
y= jLrndCTr/Srf' = 49 :r{LTPsnu(TT)Pf , 



' A(LTl3snu(JT)-, 



(10) 



lOOeV 

is much larger than one for LjjSsriuCrT ^ 
l/4(r/100ey)-'/32_j. gT.eff. includes therefore all pho- 
tons produced down to an energy that allows them to be 
upscattered to T^. For Bremsstrahlung emission, which we 
assume to be the main source of photons, the number of pho- 
tons generated diverges logarithmically at low energy, so that 
Qi.eS. may be significantly larger than the Bremsstrahlung 
generation rate of photons at T^. 

In order for the photon energy to be significantly increased 
by scattering, it must be scattered ^ nieC^ /{AT) times before 
getting re-absorbed. This sets a lower limit to the frequency 
of photons that should be included in Q^.^s., 
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where n = ln„ = 7 x 10'^«i5 cm"^. The generation rate should 
include all the photons that are generated with energy above 
hi^a and are able to be upscattered during the available time 
Lj/ilic). The Bremsstrahlung effective photon generation 
rate is thus given by 



27,eff = aetlpHeaTC 



-Aeff^eff, 



(12) 



where g^s is the Gaunt factor and Aeff ^ log[r/(/ii^n)]. As 
long as the Compton y parameter is large y > \og{T /hva), all 
photons produced above Va are upscattered to ~ 3r energies. 
Using equations (fTZt and (|9]l we obtain 



lOOa^Aeffgeff y/meC^TnipC^ 

15/8 -1/8.-1 -1 
16e, Aeff^eff- 



(13) 



(14) 
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This implies that for large shock velocities, 

/3. > 0.07«;f °(Aeffgeff)'*/'^ 

the length required to produce the downstream photon density 
is much larger than the deceleration scale. For lower shock 
velocities, thermal equilibrium is approximately maintained 
throughout the shock. 

2.4. Immediate downstream 

The photons that stop the plasma at the deceleration region 
must be generated within a distance ^ (n„/3jCT7-)~' of the de- 
celeration region in order to be able to reach it (before being 
swept downstream by the flow). Consider the following sim- 
ple model for the first few diffusion lengths downstream of 
the deceleration region. We assume that the velocity is /3rf 
downstream from some point x = xq and is much larger up- 
stream of it. As the e-folding distance of the residual velocity 



Sf3 = l3-f3d is very short, = Psnu^TT /21, this is a good a approx- 
imation. We neglect any generation of photons upstream of 
the transition, Q^^fix < xq) = 0. This is reasonable as the con- 
tribution of photons from a slab, one diffusion length wide, is 
roughly proportional to f3~^ and /? changes by a factor of 7 
along the transition region. We further assume that the gener- 
ation rate of photons downstream of the transition is roughly 
constant and equal to Q^.^tfix = xq). The latter simplification 
is justified, to an accuracy of a few tens of percents, due to 
the fact that the temperature changes by order unity across 
the available (diffusion length) distance and Q,,// oc T*^ ^ pgr 
such a flow, the density of effective photons at jc = xq is: 
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As the flow at this point is close to the downstream velocity, 
we have n^(x())Ts = p^j where T^ = r(xo) is the temperature 

the 



in the immediate downstream. Using Eqs (fT2l l and ( ITsI l 
momentum conservation equation can be written as 



^a^.AeffgeffV' 



meC^TsHdlSj =-nul3^mpC 



(16) 



where Aeff = lOAeff.i ^ \og{Ts/{hvci)\ with hva given by 
Eq. (fTTI) . Here we assumed that the temperature does not 
change much within a distance of Ljec downstream of the 
point X = X(). The Compton y parameter, relevant for the up- 
scattering of the photons that diffuse up to the deceleration 
region, 



AZ 



Pf ~ 400 



lOkeV 



PS. 



(17) 



is much larger than 1 for ^ 30/3^_[eV. We therefore as- 
sume that all photons above hvg are upscattered all the way to 

Solving Eq. ( fTSI l for /3j we find: 



A = Qa^Aeff^eff 
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Eq. (fTSI) is in good agreement with the results of (IWeaverl 
IT976h for Ts < 50 keV. For /3, > 0.07 and n„ = 10'^ cm""* - 
10^^ cm~^ the temperature (velocity) given by equation ( fTsT i 
for a given velocity (temperature) agrees with the numerical 
calculations to within a factor of 1.5 (20%) (with the gaunt 
factor calculated as in Weaver 1976). We note that the scaling 
with parameters in equation ( fTSl l can be found by wr it ing th e 
RMS equations in dimensionless form (cf Eqs. (lAl lb -( IA16l l). 

Comparing Eq. (fTSl l with the equation for the far down- 
stream temperature Td, eq. (|5]l, we see that the shock tem- 
perature is much larger than the downstream temperature for 
shock temperatures Tj > 1 keV corresponding to shock veloc- 
ities (3s 0.1 (see also Weaver 1976). 

Note, that for temperatures T >\ keV the typical photons 
in the immediate downstream have energies hv ^ 3T higher 
than the maximal energies attainable by bulk Compton accel- 
eration in non-relativistic flows, given by Eq. Indeed, 
the typical photon energies, as derived from ( fTsT l. rise much 
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faster with velocity than /jj/max oc 0^, and at T ^ 10 keV, cor- 
responding to fis ~ 0.2, are already much higher {hv ^ 3T ^ 
30 keV compared to hv^a^ ^ 4 keV). At temperatures exceed- 
ing r ^ 50 keV, where pair production becomes important, 
Eq. ( fT6] l is not applicable anymore. However, the velocity re- 
quired for bulk Comptonization to be able to accelerate to the 
typical photon energies at those temperatures, hi/ ^ 150 keV, 
approaches c. Thus, for non-relativistic and mildly relativis- 
tic shock velocities with /J^ > 0.1, bulk photon acceleration 
cannot reach energies that are considerably higher than the 
typical "thermally"-Comptonized photon energies in the im- 
mediate downstream. 

2.5. Description of the shock structure 
We can broadly divide the shock into four separate regions. 

1. Near upstream: A few diffusion lengths, (/3s(T7n„)"', 
upstream of the deceleration region. In this region, 
characterized by velocities that are close to the up- 
stream velocity, (3 w /3j, and temperatures T 3> T,,, the 
temperature changes from r„ to ^ 7^. It ends when the 
fractional velocity decrease becomes significant. 

2. Deceleration region: A (/3st77-n„)"' wide region where 
the velocity changes from /3„ to Pd and the temperature 
is roughly constant, T criT,. 

3. Immediate downstream: Roughly a diffusion length, 
{PsCTtiIuT^ downstream of the deceleration region. In 
this region, characterized by velocities close to the 
downstream velocity, (3 w f3d, and temperature T ^ T,, 
the photons that stop the incoming plasma are gener- 
ated. Upstream of this region (3 > /3d and the photon 
generation rate is negligible. Photons that are generated 
downstream of this region are not able to propagate up 
to the transition region. 

4. Intermediate downstream: The region in the down- 
stream where most of the far downstream photons are 
generated and T changes from to Td- This region 
has a width Lj given by Eq. (T3[ . much greater than 
(PsCTrnuT^ ■ Thus diffusion within this region can be ne- 
glected. The temperature profile is expected to follow 
T cx x~^. To see this, note that the photon density at a 
distance x from the shock is proportional to the integral 
of the photon generation, cx T~^l^x. Since the pho- 
ton pressure equals the downstream pressure, we have 

oc T"' and T cx x~^ (this is valid for a constant value 
of Aeffgeff and is somewhat shallower in reality). Using 
this dependence of the temperature on distance, we find 

that LTl3snu(7T ^ Aeffgeff Uj^^^ [(Aeffgeff) liTT^'^^]"' , in 
agreement with equations (fl3T l and (fTsT l. 

As an illustration, the profiles of the velocity and tempera- 
ture of a non-relativistic RMS (f3s = 0.25, n„ = 10'^ cm"-^) are 
shown in figure [1] A description of a novel method used to 
find such solutions is described in § |A] As can be seen, the 
analytic estimates are in good agreement with the numerical 
calculation. 

3. RELATIVISTIC RADIATION MEDIATED SHOCKS 

3.1. General discussion 

The analysis presented above is valid for temperatures up 
to about Ts ~ 30 keV, where a number of effects that were 




(x-xmin)(3^n B 

Fig. 1. — RMS Shock temperature and velocity profile for e = 30MeV 
(Bs = 0.25) and «„ = lO" cm"^^. The black full line is the velocity given by 
(8). The Blue full line is the temperature calculated numerically [see §|A|. 
The dashed green line is the far downstream equilibrium temperature given by 
Eq. (5)- The dashed-doted red line is the analytic estimate for the immediate 
downstream temperature. The width of the temperature equilibrium length, 
given by eq. )13t . for these parameters is L7/3„(77-n„ ^ 130. 



not taken into account begin to affect the temperature pro- 
file. The main process that should be taken into considera- 
tion is the creation of electron-positron pairs at T > 50 keV. 
Additional corrections include relativistic corrections to the 
Bremsstrahlung emission rate, inclusion of photon creation 
by double Compton scattering, Klein-Nishina (KN) correc- 
tions to the Compton cross section, reduction of the Compton 
y parameter due to the smaller optical depth, r ~ /3j^, and the 
anisotropy due to the high velocities, (3,, ^ I. 

The following important characteristic of the shock struc- 
ture is likely unchanged: the plasma decelerates due to scat- 
tering with photon and pairs that are produced in the first (3'^^ 
optical depths of the downstream [a discussion of the decel- 
er ation velocity profile, n e glecti ng KN corrections, is given 
in iLevinson & Bromberd (120081) . note that KN corrections 
and significant electron-positron production significantly af- 
fect the result (Budnik et al. 2010)]. As the downstream ve- 
locity is f3d < I /3 for all /3s, the diffusion approximation in 
the downstream is reasonable and an equation similar to Eq. 
( fTSI l. with the relativistic corrections included, can be used to 
constrain th e temper a ture in the transition region. 

Note that iWeaveii (119761) claimed that once a significant 
amount of pairs is created, the pairs will 'short out' the elec- 
trostatic field which is responsible for the deceleration of the 
protons. He claimed that the electrons and ions will no longer 
be fully coupled, and suggested that once there are many pairs 
per proton, there will be a relative drift of order between the 
leptons and protons, and that the main mechanism for stop- 
ping the protons will be Coulomb collisions with the leptons. 
We find this picture to be unrealistic in real plasmas. A rela- 
tively weak magnetic field, which will be freezed into the pair 
plasma, will suffice to deflect the protons on scales which are 
much shorter than determined by Coulomb collisions. Sev- 
eral length scales may be important here. The gyroradius of 
the protons in a magnetic field B = 10^B2G is 

2 

Rl = 7/3^ ~ 3 X IOV/3B2' cm, (19) 
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the skin-depth of the pair plasma is 



•sd ' 



1.7 X 10" 



-1/2 

'Ills cm, 



and the Compton mean free path is, 
1 



riiUT 



1.5 X 10 «/ 15 cm. 



(20) 



(21) 



As long as the Compton mean free path is much larger than 
the other two scales, it is probably safe to assume that the pro- 
tons and electrons are highly coupled on the dynamical scales. 
We thus assume in what follows that collective plasma effects 
will keep the pair plasma and the protons coupled without the 
need for an electrostatic field. 

3.2. T-j3j relation 

As long as the Compton y parameter is much larger than 
unity, the spectrum will be close to a Wein spectrum. For the 
entire parameter space considered here, the positron number 
is in equilibrium. To see this, note that the typical time it takes 
a positron to annihilate is similar to the typical time it takes a 
photon to Compton scatter, thus pairs can be annihilated (or 
generated) in Compton optical depths compared with the 
available. 

A detailed discussion of the emiss ion of pla s mas i n pair 
and Compton equilibrium is given by ISvenssoiil (11984) . who 
shows that the main emission processes are Bremsstrahlung 
and Double Compton, with a domination of Bremsstrahlung 
at temperatures above T « 60 keV. 

A precise determination of the relation between and /3„ 
is complicated in the range Tj > 30 keV due to the relativistic 
corrections to the emission rates and due to the order unity y 
parameter Such determination requires a careful analysis of 
these effects. We avoid such a detailed analysis by limiting 
our analysis to a demonstration of the validity of the follow- 
ing statements, which set stringent constraints on the relation 
between 7^ and /3d : 

1 . Temperatures 71 > 30 keV are reached for downstream 
velocities (3^ > 0.03, with « 0.03 for 7; « 30 keV 
(as shown in §|2]i; 

2. At (3^ = Q.l the shock transition temperature satisfies 
60 keV < T,; 

3. Ts < 200 keV for all the possible downstream veloci- 
ties, < 1/3; 



T he v alidity of statements 2 and 3 is demonstrated in § l3.3l and 
§ l3.4l below by comparing the photon to pair ratio that is de- 
termined by the balance of photon production and diffusion, 
eq. ( fTSl ). with the ratio that is determined by pair production 
equilibrium. The complications due to the various corrections 
are avoided since, under the conditions we consider, the fol- 
lowing hold: The number of pairs greatly exceeds the number 
of protons, n; ^ Hp, where ni = n+ + n-K, 2n_; Double Comp- 
ton emission is negligible; The Compton y parameter is large. 

3.3. At PdTuQ.l, 60 key < T, 

Assuming pair production domination and neglecting Dou- 
ble Compton emission, Eq. ( fTSl l can be written as 



n 1 

— = TaeAeff^eff,rel(r)/3j^, 

«/ 3 



(22) 



where we wrote the free-free emission in the form 

g^^eff = aeO-7-Cn^Aeff^eft;rel(r). 



(23) 



Here, t = T /(m^c^), and ^eff.iei is the total Gaunt factor [de- 
fined by Eq. ( |23]) 1 including all lepton-lepton Bremsstrahlung 
emission. For 10 < Aeff < 20 and 60keV < T < nieC^ the 
approximation 

^eff,rel ~ Aetf/2, (24) 



agrees with the results of lSvenssoiil (119841) to an accuracy of 
better than 25%. Substituting Eq. in Eq. (|22l) we find 



(25) 



We next estimate the absorption frequency. The absorption 
is dominated by e + e- Bremsstrahlung and for T > QAnigC^ 
photons require a single scattering to change their energy con- 
siderably. The absorption frequency is thus roughly equal to 
the threshold frequency i'c=b at which the Compton scattering 
rate equals the absorption scattering rate and is given by 



hl^a ~ hl'c=B = nieC 
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where «/ = 10'^n/_ig. The lepton number density «; is roughly 
given by 



n, = 
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2^2 



111 \ IflpP^C 



50keV 



(i^nxs cm ^ 



yielding 



A„=ln(-^)«ll-iln(«i5), 



(27) 



(28) 



where we used ^ 10,n;/n-y ^ 10 and T ^ 50keV (the 
result is not sensitive to the specific values adopted for these 
parameters). 
The Compton y parameter is 



y = (4r+ 16f ^)/3^2 > 40(r/50 keV)/3^-ii. (29) 



Equations ( |28] | and ( |29] l imply that for T ^ OAnigC^, the 
Compton y parameters is large enough, y > Ag, to allow all 
the photons above the absorption frequency, hva, to be comp- 
tonized (photons are upscattered by a factor of e'' in energy 
moving all photons above hi^a to the Comptonized energies 
3T = 3e^"lwa) . This implies, in turn, that photon production 
at all frequencies above hi^a contribute to the effective gener- 
ation rate and Aeff = A^. 

Let us compare now the photon-lepton ratio derived for the 
balance of photon production and diffusion, to the ratio de- 
rived from pair production equilibrium. For large chemical 
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potential, fich ^ T pair production equilibrium gives 



ni 



|eq ■ 



dxx e 



,2 _ ^2(r ') 



7^2 



(30) 



where K2 is the modified bessel function of the second kind 
of order 2. At the temperature range 60 keV < T < 90 keV 
this ratio changes by an order of magnitude: n-y/n/leqC?" = 
60 keV) « 130 while n^/n,|eq(r = 90 keV) « 13. Compai'ing 
these values to the result given by equation dZSl l. we conclude 
that at = 0.1 the temperature must be above 60 keV. 

3.4. The relativistic limit /3d = 1/3 

In the regime 200 keV < r < niec^, the pair production 
equilibrium is approximately given by 

n^/ni fi^O.Sm^c^/T. (31) 
The number of leptons is now given by (see Eq. (l27Ti). 

rii \ T^nipC^ 



111 ' 



3T 



0.6xlO'M- 



r;ni5 cm \ 



(32) 



and the absorption photon energy by 



1/2/ T V'^' 1/2 



n^'igeV, 



(33) 



yielding 



A« = ln(-^)« 15-iln(ni5)-ln(r,), (34) 

where we used g//,;/ ~ 10,n;/M^ ~ 1 and T ^ irieC^ (the re- 
sult is not sensitive to the specific values adopted for these 
parameters). 

We can rewrite equation ( |25T l for f3d ~ 1/3 and 
60 keV < r < niec^ as 



rii 



;2.5 



Aeff 

15 



(35) 



Comparing equations (l35T l and (IJTT i. we see that if 
T > 200 keV there would be too many photons generated per 
lepton. At these high temperatures, the Compton y parameter 
is large and radiative compton emission is negligible. At tem- 
peratures r > 100 keV KN corrections to the Compton cross- 
section of photons with energy ^ 3T become significant. 
The possible effect that KN corrections have is to lower the 
temperature due to the larger photon generation depth. This 
in turn would suppress the KN corrections. Hence, KN cor- 
rections cannot change the conclusion that T < 200 keV. We 
conclude that for > 0.1, 60 keV < T < 200 keV. 

4. SUPERNOVA SHOCK BREAKOUTS 

Once an expanding radiation mediated shock reaches a 
point where the residual optical depth is of order r ~ /?"', 
whether at the outer shell of a star's envelope or in the sur- 
rounding wind (in case of an optically thick wind), the radi- 
ation escapes, and the shock no longer sustains itself. In this 



section we summarize the relations between the stellar/wind 
parameters and the velocity, energy and duration of the result- 
ing shock-breakout radiation outburst. We first derive general 
relations betwee n the energy, velocity, radi us an d duration of 
the outburst in § 14.11 Next, we discuss in § l4.2l the dynamics 
of shock acceleration near the stellar surface, and its evolu- 
tion within the wind, and derive characteristic breakout veloc- 
ities. The X-ray characteristics of the breakout are d iscussed 
in § gj] based on the results of § H § |3] and § g^] Finally, 
the implications of our results to the interpretation of recent 
X-ray outbursts associated with SNe (XRO080109/SN2008D, 
XRF060218/SN2006aj), and possibly to other SN-associated 
outbursts (XRFs/GRBs), are discussed in § 14.41 In this 
section, we restrict our discussion to non-relativistic shocks 
f3s < 0.6 and non relativistic temperatures Tj < 50 keV, cor- 
responding to photon energies hi/ < 150 keV, where pair cre- 
ation can be neglected. 

4. 1 . A general relation between the breakout energy, 
velocity, and radius 

Consider a shock with velocity (3^ ^0.1 approaching the 
photosphere of a spherical mass distribution. The break-out 
will occur once the optical depth down to the photosphere is 

T - < 10. 

The energy carried by the photons in a shell of shocked ma- 
terial with a width of the order of the shock width is roughly 
given by (e.g. lWaxman et af]|2007h 



5E^ = fEdm^y/l 



(36) 



where 6m is the mass of the shell and /e is a coefficient that 
depends on the velocity and on the geometry. For a shocked 
shell in the downstream of a non-relativistic planar shock, 
/£ = 36/49. /e is of order unity for non relativistic and mildly 
relativistic shocks in spherical geometry, as long as the dis- 
tance to the photosphere is not much larger than the radius. 
The mass of a shell that has an optical depth r is roughly 
given by 



Smr^4TTR^TK-^ « 10^''R\2(K/KTy^T0.5g, 

where kt = crr/nip ^ 0.4cm^g"', r = IO'^'^tq.s 



(37) 
and 



R = 10'^7?i2 cm. Assuming that the energy of shocked plasma 
with optical depth t ^ 3to.5 is emitted, we find that there is a 
simple relation between the released energy, breakout radius 
and shock velocity at breakout, 

SE - 0.5Sm(P,cf - 4 X 10^^ R\2{k / ktT^to. 5 erg. (38) 

Once a significant fraction of the energy in the shock re- 
gion is emitted, the steady state shock solution is no longer 
applicable. Note that before the solution becomes invalid, a 
significant fraction of the energy in the shock region must be 
emitted. We thus expect that photons with energies of the or- 
der of 3Ts carrying an energy similar to the energy estimated 
in ( |38] | are inevitably emitted. 

4.2. Breakout Dynamics 

In this subsection we discuss the dynamics of shock accel- 
eration near the stellar surface, and its evolution within the 
wind, and derive characteristic breakout velocities. Many 
of the relations derived can be found in the literature (e.g. 
iMatzner & McKee|[T999l) . 

At the outer layers of stars, where the enclosed mass, M, is 
approximately constant, M = M* = const, and the composition 



is approximately uniform, the density distribution is a power- 
law in the distance from the edge of the star, 

P = P,5'\ (39) 

where 7?* is the stellar radius and 

5 = {R,-r)/r. (40) 

This relation is exact if we assume M = and a polytropic 
equation of state p = kp'^", in which case n = (7„ - 1)" 

■IGMV 



and 



Pi 



(41) 



For an efficiently convective envelope, the entropy is con- 
stant. Assuming an adiabatic index 7„ = 5 /3, we have « = 3 /2. 
For a radiative envelope, with radius independent luminosity 
L = < LEdd and mass M = M», we have = {L^/LEdd)p, 
where L-em = 4-ttGMc /k is the Eddington luminosity limit, re- 
sulting in a polytropic equation of state p cx p"*/^, or n = 3. The 
optical depth from the edge of the star is given by 



pndr oc S" 



The density at the outer edge is given by 



P=(Pi)" 



(«+ l)r 

kR^, 



«/(«+!) 



Denoting 



Pl=f, ' 



IP 4tt 



3 



(42) 



(43) 



(44) 



where M^j is the mass of the ejected envelope, we find 



p « 2.0 X lQ-'fy\M,,/MQy/\n/nT)-''%l''4!,\ 



gem 



for n = 3, and 



(45) 



p « 2.8 X lO-'YAM^>/Mef''(^/^Tr'^%T4'h cm"^ 

(46) 

forn = 3/2. 

As the shock wave approaches the stellar surface, it acceler- 
ates down the density gradient. In planar geometry, with den- 
sity vanishing at x = xo following p oc [ixo-x)/xo]" = S", the 
flow appr oaches a self-similar behavior as the shock velocity 
diverg es (iGandeTMan & Frank- Kamenetskiil [79561: ISakurai I 
Il96d) . In fliis limit 



v _5 CX S 



(47) 



where Ay = /?!«/(«+ 1), and /3i ~ 0.19 for 7 = 4/3 and n = 
3,3/2. The planar self-similar solution provides a good ap- 
proximation for the shock behavior near the surface of the 
star, where the thickness of the layer lying ahead of the 
shock is small compared to the radius, i.e. for 6-^1 (e.g. 
iMatzner & McKeelfl999h . Denoting the mass M, -M(r) of 
the shell lying ahead of radius r, at 5 < (/?, -r)//?», by 
m = M» - M(r), and noting that near the surface the optical 
depth T is r cx m, equation (l47T i may be written as 



V j cx m 



(48) 



with A,, = /?!«/(«+ 1). The kinetic energy carried by the ma- 
terial having velocity greater than some v^ scales with the ve- 
locity as 



E(> V,) cx Vj 



(49) 



where 



(5„ = 



1 + ^2A 

/3i 



for n = 3 and n = 3/2 respectively. 



^5,7 



(50) 




Fig. 2. — Shock velocity amplification, and density typical of the edge of a 
Blue Super-Giant {n=3,MseS =Mq,R = lO'^ cm ), as a fu nctio n of the optical 
depth from the edge of the star as given by Eqs. j52) and i45\ . The density is 
normalized to the value at t = I . 

The velocity of the shock at an optical depth r is larger than 
the typical bulk envelope velocity, Vgj.b. ~ y/E^j/M^. where 
Ef.j is the energy deposited in the envelope, by a factor 



V 4^«2 



(51) 



where Mses is the mass of the region in which the planar so- 
lution is applicable. For n = 3 we have Av ~ 0.14, while for 
n = 3/2 we have Ay « 0.11. The velocity amplification factor 
for n = 3 is 



ll(Mse./MQ)"-^^^/^^)"'X2'Vol'\ 



and for n = 3/2 it is 



4{M,,,/MQf'\K/KTr'R^','\ 



0.11 D-0.23 -0.11 



0.5 



(52) 



(53) 



The density and velocity profiles of the outer envelope are 
illustrated in figure |2] 

If there is an optically thick wind surrounding the star, 
breakout occurs within the wind. Assuming a wind with a 
constant mass loss rate M and velocity v,,,, p^. = M/47rv,v7?^, 
the optical depth between R and infinity is related to the mass 
m,,(< R) by r(> R) = m,,(< R)k/AttR^ (for R > /?*), and 
breakout is expected to occur at a radius 



Rh, = 



M K 
V,,. 4ttt 



'0.6x 10' 



M_5 
V,v,8 



(k/k7-)T()5 



(54) 



At early time, the velocity Vj u' of the shock driven into 
the wind at radius R is approximately given by the veloc- 
ity of the fastest shell of the ejecta. Since the final ve- 
locity reached by a fluid element shocked at v,, following 
its accelerati on due to the post-shock adiabatic expansion, 
is V/ ~ 2vj ( IMatzner & McKed [19991) . the velocity of the 
shock driven into the wind is initially 2Av,Vej. At later 
time, when the accumulated wind mass begins to decelerate 
the ejecta, the shock velocity is approximately obtained by 
equating the mass of the shocked wind, m„{< R), with the 
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mass of the part of the ejecta that was accelerated to velocity 
V > V.5,,,,,. Using eq. (|48T l we have m„(< R) w Mejiv^.w/^yeiT^" 
or v^_,^,/vej w 2[Mej/mH,(< R)]^\ This gives: 



2Av ,30 



Ma 



for BSG (n = 3) and 



Ver 



2Av.J3 ( -2^ 



1/7 



1/9 



M_5 
V„.,8 



M_5 
V>,,8 



-1/7 



-1/9 



7? 



■1/7 



■1/7 



(55) 



(56) 



forRSG (« = 3/2). 

The assumption that the shock velocity in the wind is equal 
to the velocity of the shocked ejecta is correct as long as the 
density encountered by the reverse shock in the ejecta is much 
larger than the density encountered by the forward shock in 
the wind. 

Assuming typical bulk envelope velocities of Vej,b ~ 

y/E~/M~, - 3 - 10 X 10** cm/s, or /3ej, b 0.01 - 0.03, the 
typical shock velocities at the last few optical depths of BSGs 
and WR stars may reach values of /S, > 0.2, for which tem- 
peratures exceeding T ^ 10 keV may be reached in the shock 
transition region. The shock velocities in the wind may be 
up to twice higher than the maximal shock velocity achieved 
within the star 

4.3. Breakout X-ray characteristics 

Let us consider the radiation that may be emitted during the 
break-out of a shock from the stellar surface. We consider 
the emission during a time interval R^,/c following the time 
at which the shock reached the surface (as long as the star is 
not resolved, any emission will be spread in time over this 
time scale). On similar time scales, the outer part of the star 
expands significantly, over a distance /3/7?* ^ R^,, where f3f « 
2/3, is the final velocity reached by a fluid element shocked at 
f3s (following its acceleration due to the post-shock adiabatic 
expansion, e.g. iMatzner & McKeelll999 ). On the same time 
scale, /?,/c, photons that originated at an optical depth Tesc ^ 
(5"'''^(t = 1) are capable of escaping. Using equation ( |42] | we 
have 



4-ttRI 



(57) 



This implies 



Tesc ~ lO(Mej/M0)l/^K/Kr)''''^ll^' (58) 

for BSG typical parameters, and 



15(Mej/Mo)i/5(K/Kr)''''/? 



-2/5 
13 



(59) 



for RSG parameters. For 1 > /3j > 0.1, these values are not 
much larger than the shock widths, /37', and a non negligible 
fraction, /J^' / Tesc, of the photons in these regions are likely 
to escape over ^ R^,/c. The amount of energy emitted during 
/?,/c is thus not much larger than the estimate given in eq. 

Note that the radiation observed at a given time is the sum 
of radiation emitted from different positions of the star at dif- 
ferent stages of the breakout (up to differences of /?,/c). Pho- 
tons that originated from shocked material in deeper layers, 
T ^ /3~\ will have characteristic energies which are lower 



than the photon energies at the transition region hv ^3Ts due 
to photon production and possible adiabatic cooling. There- 
fore we expect to see a spectrum which is a sum of Wein-like 
spectra with temperatures reaching the tens of keV shock tem- 
peratures Ts and starting from lower energies. 

Let us consider next breakouts that take place within an 
optically thick wind. In this case, the optical depth of the 
shocked wind material, which is compressed to a thin shell, is 
comparable to the optical depth of the wind ahead of it. Thus, 
we expect all photons contained within this region to escape 
on a time scale similar to /?/Vs, the time scale on which the op- 
tical depth changes significantly. Again, the spectrum would 
be composed of Wein-like spectra with temperatures reaching 
7;. 

Finally, the following point should be stressed. As the 
shock approaches an optical depth of 1 //3 from the surface, 
the shock structure deviates from the steady state solution due 
to the escape of energy from the system and due to the steep 
density gradient. Note, that at this stage the shock propagates 
within a density profile which varies on a length scale com- 
parable to the shock thickness. As for the effect of energy 
escape, a significant fraction of the energy in the shock tran- 
sition region must be emitted before the structure is changed 
and thus the steady state estimate for the photon energies will 
correctly describe the photons carrying most of the escaping 
energy. It is reasonable to expect that the effect of the den- 
sity gradient on the temperature will be of order unity for fast 
shocks, as the pre-shock density across a shock width changes 
only by a factor of few- by ^ 1 //? in the case of stellar surface 
[see Eqs. ( l45] ll and by ~' 1 in the case of a wind. An exact cal- 
culation of the emerging spectrum requires a time dependent 
calculation and is beyond the scope of this paper. 

4.4. Implications to recently observed SN X-ray outbursts 

It was long been suggested that a strong outburst of Ul- 
traviolet or X-ray radiation may be o bservable whe n the 
RMS reaches the edge of the star (e.g. 'Colgate! [1974 : 
1 19781: iKlein & Chevalier 1978; Ensman & Burrows! 
Matzner & McKeQil999i: iBlinnikov et al.ii200 0). These out- 
bursts carry important information about the progenitor star 
dArnettj Il977h including a direct measure of the star's ra- 
dius. During the past two years, the wide field X-ray de- 
tectors on board the Swift satellite detected luminous X-ray 
outbursts prece ding a SN explosion i n two cases, SN2006a j 
and SN2008D (iCampana et al.ll2006t ISoderberg et al.ll2008h . 
Analysis of the later optical SN emission revealed that 
both were of type Ib/c, probably produced by com- 



If^ 

1991 



pact 



( WR) progenitor sta r s (jPian et al 



2006; 



2006 t lModiaze talJ 120061; iMaedaetal.. 
20081 Malesani ~lJ 120091; iTanaka et all 1200811 

authors argue that the 



Mazzali 
Mazzali 



et aU 
et all 



While 

X-ray outbursts are produced 
(Campana et al. 2006; Waxmanet al] 



some 

by sh ock breakouts 

I2007t ISoderberg et alJ |2008), others argue that the spec 
tral properties of the X-ray bursts rule out a breakout in- 
terpretation, and imply the existence of relativistic ener- 
getic jets penetrating through the ste llar mantle/envelope 
(Soderberg et al. 2006; Fan et al. 2006; GhiselUni et"ani2007l; 
IIJ I2007t iMazzali et all l2008t iLj ,2008) (some argue that 
the break out interpretation holds for SN2008D, but not for 
SN2006aj lChevaUer & Franssonll2008h . 

The main challe nge raised f or the breakout interpre- 
tation (e.g. Mazzali et al.l 120081) was that in both cases, 
SN2006aj and SN2008D, the observed X-ray flash had a non- 
thermal spectrum extending to > 10 keV, in contrast with the 
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sub- k eV thermal temperatures expected dMatzner & McKe^ 
119991) (Note, that while thermal compone nts with low temper- 
atures or fluxes can not be ruled out, e.g. iMazzah et alj|2008t 
iModjaz et al.ll2009l such thermal components carry only a 
small fraction of the observed flux). However, the analysis 
presented here implies that fast, /3s > 0.2, breakouts may be 
expected for compact, BSG or WR progenitors, and that for 
such fast breakouts the X-ray outburst spectrum may naturally 
extend to lO's or lOO's of keV. 

Let us discuss in some more detail the possible interpre- 
tation of the X-ray outbursts associated with SN2006aj and 
SN2008D as breakouts. In the case of SN2008D, both the 
X-ray outburst and early UV/O emission tha t followed it are 
consistent with a breakout interpretation (see lSoderberg et aTl 
l2008i for detailed discussion). The X-ray burst energy and du- 
ration are consistent with the relation of eq.|38] and the UV/O 
emission is consistent with that expected from the post-shock 
expansion of the envelope. Moreover, the non-thermal X-ray 
and radio emission that follow the X-ray outburst are consis- 
tent with those expected from a shock driven into a wind by 
the fast shell that was accelerated by the radiation mediated 
shock, for wind density and shell velocity, v/c ^ 1/4, and 
energy, E ^ 10'*^ erg, which are inferred from the X-ray out- 
burst. The fast velocity inferred for the late radio and X-ray 
emission is consistent with that required to account for the 
non-thermal spectrum. 

The case of SN2006aj is more complicated. The X-ray 
burst energy and temperature are consistent with a mildly rel- 
ativistic, v/c ^ 0.8, shock breakout from a wind surrounding 
the star, the UV/O emission is broadly consistent with that 
expected from the expanding envelope, and the non-thermal 
X-ray emission that follows the X-ray outburst is consistent 
with that expected from the shock driven into the wind by 
the fas t shell that wa s accelerated by the radiation mediated 
shock dCampana et al. 2006; Waxman et al. 2007). However, 
the duration of the X-ray outburst is larger than expected, 
the UV/O emission is not as well fit by the model as in the 
case of 2008D, and the non-thermal radio emission is higher 
than expected from the wind-shell interaction. Several au- 
thors (Campana et al. 2006; Waxman et al. 2007) have sug- 
gested that the deviations from the simple breakout model are 
due to a highly non-spherical explosion. Some support for 
the non spherical nature of the explosion was la t er ob tained 
the polarization measurements ofGorosabel et al.'('2006) (see, 
however, Mazzali et al. 2007). Other authors (Soderbera et al. 
120061: GhiselHni et al. 2003 IFan et al.ll2006t ILi et al.ri2007n 
have argued that a relativistic jet is required to account for 
the observations. We believe that an analysis of the modifi- 
cations to the simple spherical model, introduced by a highly 
non-spherical breakout, is required to make progress towards 
discriminating between the two scenarios. 

There is an additional important point that should be clar- 
ified in the context of SN2006aj. Under the shock breakout 
hypothesis, the energy and velocity of the accelerated shell 
that is responsible for the X-ray emission are v/c ^ 0.8 and 
E ^ lO'*'''^ erg. This is similar to the parameters of the fast 
expanding shell inferred to be ejected by SN1998bw (asso- 
ciated with GRB080425), based on long term radio and X- 
ray emission, which are interp reted a s due to inter action with 
a low density wind (Kulkarn i et al. |[l998; Wax man & Loebl 
[T999t lLi&Chevaliei1il999i: iWaxmanl l2004a.ibi) . Long term 
radio observations strongly disfavor the existence of an en- 
ergetic, 10^' erg, relativistic jet associated with SN1998bw 



dSoderberg et al.1 120041) . The similarity of SN1998bw & 
SN2006aj may therefore suggest that such a jet is not present 
also in the case of SN2006aj. However, the large energy de- 
posited by the explosion in a mildly relativistic shell, v/c ^ 
0.8, is a challenge in itself: The acceleration of the super- 
nova shock near the edge of the star is typically expected 
to de posit only ^ 10^^ ergs in such fast shells (iTan et al.l 
l200l . as can be inferred, e.g, from eq. ( |49] l, which implies 
£(> v,)/£ej. (v./vej.,b)-^^ 

5. DISCUSSION 

We presented a simple analytic model for non-relativistic 
(§ |2]l and relativistic (§ O radiation mediated shocks. At 
shock velocities ft = v,/c > 0.1(«/10'^ cm'^)'/^" [Eq. (Ell], 
where n is the upstream density, the plasma is far from ther- 
mal equilibrium within the transition region, where most of 
the deceleration takes place, since the plasma does not have 
enough time to generate the downstream black-body pho- 
ton density n^.-y ^ obbTJ /3. The electrons (and positrons) 
are heated in this region to temperatures significantly ex- 
ceeding the far downstream temperature, ^ T4, and the 
photons are in Compton equilibrium with the electrons (and 
positrons). The transition temperature, 7^, is independent of 
the upstream density: It is given by /3s ~ 0.2(7^/10 keV)'/** 
for /3s < 0.2 [eq. At velocities > 0.6 (^d > 0.1), 

where the plasma is dominated by electron-positron pairs in 
pair production equilibrium, the temperature is constrained 
to the range 60 keV < 71 < 200 keV. The thermalization 
length, the distance behind the transition region over which 
the plasma thermalizes and reaches the downstream temper- 
ature Td, is ^ {Ts/Tctyl^ times larger than the transition (de- 
celeration) width, ^1/ /3s<JTn. Our simple model estimates 
are in agreement with the results of exact numerical relativis- 
tic calculations, whi ch will be presented in a follow-up paper 
jBudnik et al. II2OIOI) . 

We showed in section § |4] that radiation mediated shocks 
breaking out of the stellar envelopes of BSGs and Wolf- 
rayet stars are likely to reach velocities /3s > 0.2. We thus 
expect that for reasonable stellar parameters, the spectrum 
emitted during supernovae shock breakouts from BSGs and 
Wolf-Rayet stars may include a hard component with pho- 
ton energies reaching tens or even hundreds of keV. This 
implies that core-collapse SNe produced by BSGs/WR stars 
may be searched for by using hard X-ray (wide field) detec- 
tors. The detection rate of such events is significantly different 
than is inferred assuming a thermal emission spectrum (e.g. 
ICalzavara & Matzneill2004 ). due to both the modification of 
the intrinsic spectrum and to the reduced absorption of high 
energy photons. A quantitative estimate of the increased de- 
tection rate is beyond the scope of this paper. As the escaping 
photon energies are highly sensitive to the shock velocity, the 
rate is sensitive to the unknown high-end part of the distribu- 
tion of the shock velocities during breakouts. 

We have argued in § 14.41 that the X-ray outburst 
XRO080109 associated with SN2008D is most likely due to 
a fast breakout: The energy, ^ lO'*^ erg, duration, ~ 30 s, and 
fast ejecta velocity ~ 1 /4, are all consistent with expected 
breakout parameters, and the hard X-ray spectrum is a natural 
consequence of the high velocity, which is independently in- 
ferred also from later X-ray and radio observations. Our anal- 
ysis shows that the spectrum of the X-ray flash XRF060218, 
associated with SN2006aj, might also be explained as a fast 
breakout. However, the breakout interpretation of this event 
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is challenged by the long duration of the X-ray flash, and by demonstrates that fast breakouts may indeed produce non- 

the high energy, E ^ lO'*''-^ erg, deposited in this explosion in thermal spectra extending to lOO's of keV, supports the viabil- 

mildly relativistic, v/c 0.8, ejecta (see detailed discussion ity of the breakout interpretation of low-luminosity Gamma- 

in §|4^41i. ray bursts/X-ray flushes associated with SNe. A major chal- 



Wan g et alJ (|2007) have suggested, based on the break- lenge for such a scenario is constituted by the requirement of 
out interpretation of the X-ray outbursts of SN2006aj and large energy deposition, E ^ 10"^^^ erg, in the fastest, mildly 
SN2008D, that all the low-luminosity Gamma-ray bursts/X- relativistic (v/c ^ 0.8), part of the expanding ejecta. 
ray flushes associated with SNe, which have smooth light 
curves and spectra not extending beyond few 100 keV (like 
those associated with SN1998bw, SN20031w, SN2006aj), are 

due to shock breakouts, and do not require the existence of This research was partially supported by Minerva, ISF and 

energetic highly relativistic jets. The present analysis, which AEC grants. 

APPENDIX 

SOLVING THE NON-RELATIVISTIC RMS EQUATIONS 

We describe below our numerical solution of the NR RMS equations. The derivation given below offers several advantages 
compared to the derivation given by Weaver ( 1976): The equations are written in a dimensionless form, making the dependence 
of shock structure on parameters (in particular, on plasma density) explicit, and the diffusion equation is solved analytically, 
making the numerical integration procedure much simpler. 

We use the following approximations: 

• We neglect the pressure of the electrons and protons; 

• We assume that the radiation and the electrons are in Compton- equilibrium with equal temperatures = Te\ 

• We neglect the corrections to the Compton scattering cross section and assume it to be equal to the Thompson cross section; 

• We assume that the process generating the photons is Bremsstrahlung emission, as given by Eq. ( fT2] l; 

• We assume that the photon distribution is described by the diffusion equation. 

Under these assumptions, the equations determining the steady state shock profile are conservation of proton flux, 

= «„/?,, (Al) 



conservation of momentum flux, 
and conservation of energy flux. 



= riuPsilSs - l3)mpC^ , (A2) 



(e^ + p^)(3c-- — i = HuPsmp^- — c . (A3) 

ivipOT ax 1 

These equations, together with the relation 

e^ = 3p^, (A4) 

form a closed set of equations for the velocit y p artic le number density. Hp = n^, and photon pressure, p^, and energy, e^. 
e^, pj and n may be eliminated from eqs. ( IAlb - (IA4l l to obtain an equation for the velocity, 

~F ~ 7^ ' ^ 

ax 6^ 

where we int roduced the dimensionless variables x = 3(TTn,if3,x and f3 = fi/ P^. Equation ( IA5l l can be solved analytically (e.g. 
lWeaved[T976h . 



1 

7^" 



(1-/3)' 



(A6) 



_(7/3-l)_ 

The photon pressure is related to the temperature through 

p^ = n^T, (A7) 
where is the photon number density. T, or n^, are determined by the photon diffusion eq.. 



dx 

where 



= Qj,.s, (A8) 



c dn.^ ( 1 drir^ , 

j cfin^- —L=cfi[n^-—^—L] (A9) 

iripOj dx \ iriuPsC^T dx ' 



11 

is the effective current of photons and the generation of photons is given by Eq. ( fT2] i. 



G7,eff = aenpHfO-rC Y Aeffgeff/abs- (AlO) 

Here , /abs = l-n~J (obbT^/^) is approximately the suppression factor due to absorption of the photons jWeaveJ 19761) . Equations 
( |A8l )-( |A7l ) determine T (and n^) for given 13 and p-^. 
Defining dimensionless variables T and riy that satisfy 



ifif 



afT, (All) 



and 



equation (IA8l l may be written as 



where 



and 



Using Eq. ( IA2I ) we have in addition 



-2 

a^n„n-y, (A12) 



^ = e., (A13) 

ax 

d 

;^ = /3(n^-— «^), (A14) 

e = Aeff^eff/abs^r'7'"'/'- (A15) 

o 



f=^. (A16) 



The solution for T and is nearly independent of /3, n^, T, since these parameters appear in the equations only through 
the product Aeff geff/a bs. which is weakly dependent on parameters. In particular, equation (fTsT l satisfies the scaling given by 
Eqs. dATTT l and ( IXl2] l. 

It is straight forward to show that the Green function, G(i,i()), for the diffusion equations (IA13b and (IA14b . defined by the 
relation 



n^(f)= / G{x,XQ)Q{xQ)dxo, (A17) 

is given by 



X<Xq 



, Afr(io) 

G(fo,i)=<( , (A18) 

' X> Xq 



where, 



-(i'-.f)_ 



1 



Mx) /v /3(i') 



di'. (A19) 



The effective velocity /3eff, calculated for the velocity profile given by Eq. ( IA6I ). is shown in figure[3] 
A numerical solution to equations ( IA13l l-( lAT6b can be obtained by preforming iterations in which is calculated from T 

using Eq. ( IA15b . then is calculated using Eq. (IA17b . and then a new value of T is obtained using Eq. (IA16b . The solution 
presented in figure[T]was obtained for an interval i,nin < i < imax, ignoring photon contributions from outside of the interval, with 
-fmin = ln(me/OTp)-ln(6)/7 chosen as the point where 1 -/3 = nig /nip, and Xmax = 400. The solution is largely unsensitive to the 
choice of the interval boundary. 
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